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Recent experiments measuring the electrical conductivity of DNA molecules highlight the need for 
a theoretical model of ion transport along a charged surface. Here we present a simple theory based 
on the idea of unbinding of ion pairs. The strong humidity dependence of conductivity is explained 
by the decrease in the electrostatic self-energy of a separated pair when a layer of water (with 
high dielectric constant) is adsorbed to the surface. We compare our prediction for conductivity to 
experiment, and discuss the limits of its applicability. 



I. INTRODUCTION 

For materials where electrical conduction occurs at the 
surface, the adsorption of water molecules can have a 
very strong effect on the conductivity. A variety of ex- 
periments over the last four decades have examined the 
effect of water adsorption on surface conductivity for ma- 
terials like quartz [l| , silica gel [2| , polymer films 0] , lipid 
membranes [3] , and various ceramics p, 0] ■ In these ma- 
terials the conductivity changes by as much as six orders 
of magnitude as the atmospheric humidity is varied, in- 
creasing with humidity in most cases by a simple expo- 
nential. This dependence has been exploited to design 
humidity detectors from ceramic materials [Si]. 

In recent years, studies examining the electrical prop- 
erties of double-helical DNA and DNA bundles have 
found a similar dependence of conductivity on relative 
humidity [3, H, HI. Experimental data shows that the 
conductivity of DNA increases by a simple exponential 
over nearly the entire range of relative humidity (10% to 
100%), and spans roughly six orders of magnitude. Ex- 
periments on flat DNA films show a similar exponential 
dependence [13, Ull, [13, While the mechanism for 
electrical conduction in DNA remains controversial, it is 
widely believed to occur through the movement of ions 
along the outer surface of the double- helix [13, [3, ■ 

The similarity between conductivity measurements for 
DNA and for other surface-conducting materials suggests 
that there may be some universal law explaining the ex- 
ponential dependence of conductivity on relative humid- 
ity. As of yet, there is no generally accepted theory of 
electrical conduction on a wetting surface. In this paper 
we propose such a theory, based on the idea of unbinding 
of ion pairs. 



II. ADSORPTION OF THE WATER LAYER 

In this paper we examine the case of a uniform, 
completely- wetting surface (contact angle 9 — 0). On 
such a surface we can consider the adsorbed water to 
form a complete, flat layer rather than a set of isolated 
droplets. The case of small but finite contact angle 
< 6* ^ 1 (incomplete wetting) is discussed briefly at 
the end of this section. 



There are a number of theoretical treatments which 
yield estimates for the water layer thickness as a func- 
tion of vapor pressure p. The most cited of these is the 
so-called BET formula [l^, which has been widely em- 
ployed to describe the adsorption of the first monolayer 
of a gas onto a substrate. Despite its practicality, the 
theoretical value of the BET formula for thicker films 
remains questionable. Its description of the water sur- 
face as a "molecular pile" ignores the effects of surface 
tension, and it is known in many cases to predict sig- 
nificantly greater adsorption at high humidities than is 
actually observed [l7l[i8|. 

The most general way of calculating the water layer 
thickness is by defining the chemical potential fi — fJ-' + li-i 
of a water molecule within the layer, where /i/ is the 
chemical potential of bulk water and /i' is the component 
of chemical potential resulting from interactions with the 
solid adsorbent. In general /i' varies with the distance z 
from the solid surface such that /i'(z — > oo) = 0. The 
equilibrium water layer thickness is the distance d at 
which /i'(z = d) is equal to the free energy of conden- 
sation for a single water molecule figj: 

Fads^kBTln{p/po). (1) 

Here, po is the saturated vapor pressure and ksT is the 
thermal energy. 

The simplest way to estimate /i' is to assume that it 
is the result of Van der Waals forces between the water 
molecule and the adsorbent. If we assume the adsorbent 
to occupy the half-space z < 0, then we can estimate the 
chemical potential /i'(z) of a water molecule at z > to 
be [2CS] 



Here, is taken to mean the volume occupied by a 
single water molecule in the liquid phase {w w 3.1 A) 
and H is the (negative) Hamaker constant for the water- 
adsorbent interaction JlJ. For water-DNA interaction, 
the Hamaker constant can be estimated as H ^ —0.4 eV 
[Hil]. The Van der Waals energy of Eq. ^ is consid- 
ered a source of "disjoining pressure" , meaning that it 
works to increase the thickness of the water layer in op- 
position to the condensation free energy of Eq. ([T]) . The 
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resulting water layer thickness satisfies ^'(z = d) = Fads, 
or 



H 



6T:kBTln{p/po) 



1/3 



(3) 



a result first derived by Frenkel [2J|. For some mate- 
rials this prediction is in good agreement with adsorp- 
tion experiments [13, [l8| , while for other materials there 
are important "structural" contributions to the disjoin- 
ing pressure associated with crystalline-like ordering of 
water molecules [1^, [2^. These "structural" forces are 
highly dependent on the molecular nature of the sub- 
strate surface, and experiments probing their significance 
for materials like quartz have produced contradictory re- 
sults [l^l ■ For the remainder of this paper we will assume 
that structural forces are unimportant. 

In the case of an incompletely wetting surface with 
contact angle < 9 1, the contribution of surface 
tension to the chemical potential /i' must be considered. 
The corresponding term ~ 9^aw^ should be added to 
the right hand side of Eq. This produces a smaller 
thickness d at all humidities, and allows for finite d at 
100% humidity [H. 



III. CHARGE UNBINDING AND 
CONDUCTIVITY 

We suppose that the surface contains some two- 
dimensional concentration N of fixed charges, which for 
the sake of discussion we take to be negative (these may 
be, for example, the phosphate groups of DNA). When 
the surface is dry, these negative charges tend to be neu- 
tralized by positive ions that are tightly-bound to the 
negatives by electrostatic attraction. Depending on the 
method of surface preparation, the positive charges may 
be either protons H+ or some metal ion like Na"*" . If the 
surface was obtained by drying from distilled water, it 
will most likely be covered by 11+ ions. If the surface 
was exposed to salty water, then it will attract positive 
salt ions in order to maintain neutrality. We assume that 
any excess salt has been rinsed out. 

In order for electrical conduction to occur, some of 
these positive-negative pairs must unbind. In the absence 
of adsorbed water, this unbinding process requires a huge 
activation energy e^/ldneoa « 3.5 eV. Here, a ^ 1 A is 
the ion radius, eo is the vacuum permittivity, and the two 
charge species are taken to have charge ±e. 

The adsorption of a water layer surrounding a bound 
ion pair, however, lowers the activation energy by provid- 
ing an atmosphere of high dielectric constant. Electric 
field lines from a single unpaired ion remain preferren- 
tially within the water layer, and as a result the electro- 
static self-energy decreases. In this section we estimate 
the activation energy associated with the unbinding of 
an ion pair immersed in a water layer of thickness d, and 
from there we calculate the density of free charges and 
the conductivity. 



The Coulomb self-energy of a charge confined within 
a "slab" geometry has been solved exactly in Refs. [2l, 
[29| . Here we present a heuristic derivation to elucidate 
the structure of the electric field. We imagine that a 
single unpaired charge produces electric field lines that 
remain inside the water layer, spreading out radially in 
two dimensions, until some distance L when they exit 
the layer and behave three-dimensionally. This is shown 
schematically in Fig. [TJ Under these assumptions it is 
easy to assemble a piecewise description of the electric 
field using Gauss's Law: 



e/2TTeQKrd, 
e/47reor2. 



d/2<r < L 
r > L 



(4) 



Here, k « 80 is the dielectric constant of water, and the 
charge is taken to reside in the center of the water layer 
(r = 0). It is straightforward to calculate the Coulomb 
self- energy U of the charge: 



U 



E{rfdV 



2 2 

\n{2L/d) ^ 



A-Ke^nd 



(5) 



The distance L at which electric field lines leave the 
film is that which minimizes the self-energy: 



mj_ 
'Wl 



0. 



(6) 



Applying this condition to Eq. ([5]) leads to the conclusion 
L — Kd/2. This length plays the role of the screening ra- 
dius, truncating the two-dimensional logarithmic poten- 
tial at the distance r — Kd/2. As a result, the self-energy 
of a free charge in the water film can be estimated as 



U 



Aneond 



(1 + lnK). 



(7) 




FIG. 1: Schematic depiction of electric field lines in the vicin- 
ity of an unpaired charge within a water layer (hatched area) 
that has condensed on top of a substrate (solid area). Field 
lines remain preferrentially within the water layer, exiting 
only at a distance L ~ hid/2. 

A more accurate estimate for U can be obtained from 
the asymptotic expression for the interaction energy V{r) 
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of two charges ±e within the water film and separated a 
distance r with d -^r <^ nd [28| : 



V{r) 



2TTeoKd 



[—7 + In (Kd/r)] . 



(8) 



Here, 7 « 0.577 is the Euler constant. The self-energy 
of each unpaired ion can be estimated a,s U = {V{oo) — 
V{d))/2, or 



U 



ATTCond 



(-7 + InK). 



(9) 



Further calculations will use this expression for U. 

In order to calculate the density of free charges at a 
given layer thickness d, we define the free energy of the 
solution when some number n of the possible N ion pairs 
in a unit area are unbound. We use as the zero of free 
energy the case where none of the ion pairs on the surface, 
of total area A, are unbound. Then, if each unbound 



positive ion occupies an area ^ 
unit area can be estimated as 

F = 2nU - nkeT ln{l /tiw^) 



^ , the free energy F per 



■In 



NA 
nA 



(10) 



The first term of Eq. ([TO]) refers to the electrostatic 
self-energy of all unpaired charges (and assumes that 
n ^ eoKfcsT/e^c?, so that there is no screening of the 
self-energy U [3^). The second term corresponds to the 
positional entropy of the free positive charges, and the 
third term represents the "hole entropy" of the unpaired, 
fixed negative charges. Here (^) is used to represent the 
binomial coefficient "X choose Y" . 

We can use Stirling's approximation to simplify Eq. 
(jlOp and then impose the equilibrium condition dF/dn = 
in order to derive a relation between n and d. We arrive 
at 



N -n 



(11) 



where Ib — /ATreonkBT is the Bjerrum length; Ib ~ 7.0 
A at room temperature in bulk water. When n <^ N, Eq. 
(fTTj) reduces to 



(12) 



The two-dimensional conductivity a — j /E, where E 
is the strength of the applied electric field and j is the 
current per unit width in the direction of the applied 
field, can be determined from the density n of available 
charge carriers as 



nefj,. 



(13) 



Here, /i is the electrical mobility of the charge carriers. 
Substituting Eq. ^ into Eq. ^ gives 



e/i 



'-{-j+lnK) 



(14) 



The dependence of the water layer thickness d on hu- 
midity may be taken as an empirical relation for a given 
surface, or it may be assumed to result from Van der 
Waals forces as in Eq. ([3]). In the latter case, 



exp 



(15) 



^ (67rfcBrin(p/po)/i/)'^' (-7 + In 



If K is considered to be a constant, then the dependence 
of conductivity on humidity implied in Eq. (|16p can be 
expressed in simplified form as 



log cr 



/3(log(po/p))'/', 



(16) 



where a and (3 are positive constants. 

In general, however, the effective dielectric constant of 
the water layer depends on its thickness. The value of 
K for a water film is known to deca.y exponentially with 
film thickness to its bulk value [SlL S2l : 



K~ 80(1-6-'*/^), 



(17) 



where A is some decay length. This has the effect of 
increasing the Bjerrum length Ib oc 1/k at low humidities 
relative to its bulk value. The behavior of the dielectric 
constant in the immediate vicinity of a DNA molecule is 
not well-known, so here we use the estimate A = 3 A. 

The mobility fi in general depends on the nature of 
the charge carriers in the system. As an example, we 
can make an estimate for proton hopping between water 
molecules (the "Grotthus mechanism"). Experimental 
data s ugg ests a hopping time of about 1.5 ps in bulk 
water [33|. Using a hopping length of 2.5 A — the hy- 
drogen bond length between water and HaO"*" — gives a 
2D diffusion constant 1 x 10~^ cm^/s and a correspond- 
ing mobility 4 x 10""^ cm^/Vs. If the total density of 
fixed charges is = 10^^ A~^, then we can estimate a 
maximum conductivity of about 10~'^fl~^. 

Fig. [5] shows a typical case for the thickness d of the 
water layer, as given by Eq. ([3]), and the conductivity as 
functions of humidity. We have used changing a dielectric 
constant k as in Eq. (jl7p and the estimated Hamaker 
constant H = —0.4 eV. 



IV. FORMATION OF A WATER CAP ON A 
FREE ION 

To this point we have assumed that the water layer is 
flat, with a constant thickness d determined solely by 
the competition between the free energy of condensa- 
tion and the attraction of water molecules to the sur- 
face (Eqs. (HI) and In other words, we have treated 
the Coulomb energy as a small perturbation. This is 
a good approximation for relatively thick films, where 
the Coulomb energy is small. At humidities lower than 
about 40%, however, we are aparently dealing with water 
layers of single- molecule thickness d ~ w. Here the large 
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FIG. 2: (a) The water layer thickness as a function of hu- 
midity, from Eq. ((3|. The dash-dotted line shows ?s/2, which 
decreases as the water layer becomes thicker and k approaches 
its bulk value. 

(b) Surface conductivity as a function of humidity (solid line) , 
from Eq. (|16|) . The maximum conductivity is 10~^S1~^, as es- 
timated in section III. The dotted line is a linear expansion 
around 50% humidity. 

Coulomb energy could force the water film to swell above 
a free ion, forming an additional cap of water as depicted 
schematically in Fig. [31 

In this way the Coulomb field of free charges con- 
tributes to the disjoining pressure on the water film. The 
creation of a water cap is resisted by the combination 
of free energies in Eqs. H]) and ([2]), integrated over the 
number of water molecules that comprise the cap. The 
surface tension energy associated with the "bulge" in the 
water surface must also be taken into account. We as- 
sume that the cap has some characteristic radius R and 
height h above the fiat layer. The free energy AFcap as- 
sociated with the formation of the cap depends on R and 
/i, and can be written most generally as 

AF,„p = AC/ b.Fvdw + ^Fads + AF,„^/ , (18) 

where A[/ is the (negative) change in electrostatic energy 
associated with the formation of the cap, ^Fydw is the 
(negative) Van der Waals energy resulting from adsorp- 
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FIG. 3: Schematic illustration of a water cap around a free 
ion, formed by the adsorption of additional water molecules 
to the region of intense electric field surrounding the ion. The 
radius R and height h of the cap can be optimized to deter- 
mine the effect of the cap formation on the conductivity. 

tion of the water molecules comprising the cap, l^Fads is 
the (positive) condensation free energy associated with 
their adsorption, and AF^urf is the (positive) additional 
surface energy. Optimizing AFcap with respect to R and 
h determines the equilbrium cap size and its contribution 
to the free energy of an unbound charge. For the exam- 
ple shown in Fig. [2] at 20% humidity, we find that the 
equilibrium AF^ap ~ —OAkBT. As a consequence, the 
conductivity is increased by only about 50%, which in 
logarithmic scale would be a barely-noticeable correction 
to Fig. [3 



V. DISCUSSION 

The prediction for conductivity of Fig. ^ is in rea- 
sonable agreement with experiments P, H, U 0, H, H, 0, 
H, 1^, [13, EH, [13, Ii3| : it reproduces the exponential de- 
pendence of conductivity on humidity, particularly in the 
well-studied range of humidities p/po > 0.3, and suggests 
that the humidity can span six orders of magnitude for a 
realistic value of the Hamaker constant H. However, for 
DNA our estimates yield a conductivity that is smaller 
than observed values by one or two orders of magnitude. 
This may be the result of an unrealistically low estimate 
of the mobility based on data for bulk water. Interac- 
tions between conducting ions and the negative surface 
charges may improve the mobility. 

It should be emphasized that Eqs. and are 
only valid when free ions are exponentially rare, i.e. 
U ^ ksT. When the water layer becomes sufficiently 
thick, the self-energy of a free charge becomes compa- 
rable to ksT and the system undergoes a Kosterlitz- 
Thouless (KT) unbinding transition at d = Ib/2 [soj . 
Strictly speaking, the formulas derived in this work are 
only valid below the transition, d < and imme- 

diately after it at 2d/lB — 1 < (InK)^^ (see Fig. 2 in 
Ref. [30|). For larger d, we expect a rapid acceleration 
in the density of free charges n. For d > Ib, n becomes 
saturated, so that the conductivity no longer depends 
strongly on humidity. In the example of Fig. [51 we expect 
the KT transition to occur at about 80% humidity, and 
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the corresponding steep rise of n should become appar- 
ent at about 85%. We are not aware of any experimental 
evidence which demonstrates this behavior. 

For conduction through an assembly of DNA 
molecules, we expect that the conductivity should de- 
pend in general on the geometry of the assembly. In 
the case of a single double-helical DNA molecule with 
an adsorbed water layer, the electric field from an un- 
paired charge behaves one-dimensionally at distances 
larg er than half the circumference of the double-helix 
[34 [35| . As a result the electrostatic self-energy is 
larger than our prediction and the conductivity should 
be smaller. For a bundle of parallel double-helices, if wa- 
ter is absorbed within the bundle then the electric field of 
a free charge exhibits three-dimensional behavior at dis- 
tances larger than the interhelical spacing. The result is 
a smaller self-energy and a larger resulting conductivity. 
Surprisingly, experiments on individual DNA molecules 
DNA films [10, 11, 12, 13], and DNA bundles [7, 8] all 
show a similar exponential dependence of conductivity on 
humidity. The lack of a geometry-dependent difference 
remains a puzzle. 

The temperature dependence of DNA conductivity was 
studied only by the authors of Ref. [7| at about 50% 
humidity. They obtained an activation-like dependence 
a oc e~^°/'^^"^ with Ea « 0.5 eV. We cannot account 
for this behavior because of the temperature dependence 
of the dielectric constant, which in bulk water has the 
dependency k oc T~^'^* [36[. Further experimental work 
examining the temperature dependence of DNA conduc- 
tivity could help improve our understanding. 



Finally, this paper has considered only the case where 
fixed charges on the conducting surface are of a single 
sign, e.g. fixed negative charges with positive counteri- 
ons. One may well ask how much conductivity should be 
expected in the case of a charge-neutral surface, where 
both positive and negative fixed charges are present on 
the surface with their respective counterions. If the den- 
sity of fixed charges N on the surface is such that fixed 
charges are very sparse, then there is no interaction be- 
tween fixed charges and we expect that the predictions 
of this paper should remain valid. Each fixed charge may 
lose its counterion via the unbinding process described in 
this paper, and free ions contribute equally to the conduc- 
tivity regardless of sign. If TV is large, however, then fixed 
charges of opposite sign may neutralize each other and 
allow their respective counterions to unbind with very lit- 
tle energy cost. If this is the case, then counterions may 
be rinsed out during the surface preparation, resulting 
in a depletion of the conductivity. A small amount of 
conduction may still occur even in the absence of counte- 
rions through the ionization of water molecules, but the 
characterization of this mechanism is beyond the scope 
of the present work. 
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